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ABSTRACT SUBSCRIPTS 


Experimental data from jet-engine tests have indi- 
cated that unsteady blade row interactions and sepa- 
ration can have a significant impact on the efficiency 
of low-pressure turbine stages. Measured turbine ef- 
ficiencies at takeoff can be as much as two points 
higher than those at cruise conditions. Several recent 
studies have revealed that Reynolds number effects 
may contribute to the lower efficiencies at cruise con- 
ditions. In the current study numerical experiments 
have been performed to study the models available 
for low Reynolds number flows, and to quantify the 
Reynolds number dependence of low-pressure turbine 
cascades and stages. The predicted aerodynamic re- 
sults exhibit good agreement with design data. 

NOMENCLATURE 

M Mach number 

P Pressure 

P t Total Pressure 

Re Reynolds number (axial chord, inlet vel.) 

S Entropy, arc-length distance 

Uoo Free stream velocity 

X Axial-direction 

a Absolute reference frame flow angle 
P Relative reference frame flow angle 
r) Efficiency 

X Pressure gradient parameter 

v Kinematic viscosity 

p Density 

a Intermittency function 

r Turbulence intensity 

0 Momentum thickness 
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n Wall normal direction 
tt Total-to-total quantity 
$ Momentum thickness 

1 Vane inlet 

2 Vane exit/rotor inlet 

3 Rotor exit 

SUPERSCRIPTS 

* Reference quantity 

Relative reference frame flow quantity 

INTRODUCTION 

Experimental data from jet-engine tests have in- 
dicated that unsteady blade row (wake) interactions 
and separation can have a significant impact on the 
efficiency of turbine stages. The effects of these in- 
teractions can be intensified in low-pressure turbine 
stages because of the low Reynolds number operating 
environment. Measured turbine efficiencies at takeoff 
can be as much as two points higher than those at 
cruise conditions [1]. Thus, during the last decade a 
significant amount of effort has been put into deter- 
mining the effects of transition and turbulence on the 
performance of low pressure turbine stages. Experi- 
mental investigations have been performed, for exam- 
ple, by Hodson et al. [2, 3, 4, 5], Halstead et al. [6, 7], 
Qiu et al. [8], Sohn et al. [9] and Boyle et al. [10]. 
These investigations have helped identify/clarify the 
roles that factors such as the Reynolds number, free 
stream turbulence intensity, pressure gradient and 
curvature have in the generation of losses. In par- 
ticular, it has been determined that [1]: 

• At low to moderate Reynolds numbers there is 
a laminar region extending some distance from 
the leading edge. 

• The boundary layer may separate, particularly 
on the suction surface of the blade. Separation 
may occur in the form of a closed bubble, or 
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as massive separation with no reattachment (re- 
sulting in large losses). The pressure surface may 
have cove separation, and small separation bub- 
bles may exist near the leading edge. The separa- 
tion bubbles often originate in transitional flow, 
while reattachment is usually in turbulent flow. 

• The interaction of incoming wakes with the 
boundary layer often creates a convected tran- 
sitional or turbulent patch, which is trailed by a 
“calmed” region. The calmed region is a relax- 
ation region between the patch and the laminar 
boundary layer. 

In parallel to the experimental investigations, there 
have been significant analytical efforts to improve the 
modeling of transition. Examples of such efforts in- 
clude the works of Mayle [11], Reshotko [12] and 
Gostelow ei al. [13, 14]. These newer models show 
promise of providing accurate transition predictions 
over a wide rangle of flow conditions [15], although 
they have yet to be implemented into the numerical 
flow analyses used by the turbine design community. 
Some recent computational investigations of interest 
include the works of Chernobrovkin and Lakshmi- 
narayana [16], Kang and Lakshminarayana [17] and 
Huang and Xiong [18]. 

The focus of the current effort is to use a 
viscous, unsteady quasi-three-dimensional Navier- 
Stokes analysis to study the models available for sim- 
ulating low Reynolds number flows. The effects of 
Reynolds number variations on performance of low- 
pressure turbine blade rows have also been investi- 
gated. Both cascade and stage simulations have been 
performed. The numerical results have been com- 
pared with the available design data. 

ALGORITHM 

The current work is based on an extension of an 
approach developed by Rai [19]. The approach is re- 
viewed in brief here. The flow field is divided into two 
types of zones. O-type grids are used to resolve the 
flowfield near the airfoils. The O-grids are overlaid on 
H-grids which are used to resolve the flow field in the 
passages between airfoils. The H-grids are allowed to 
slip relative to one another to simulate the relative 
motion between rotors and stators. The thin-layer or 
full Navier-Stokes equations are solved on both the 
O- and H-grids. The governing equations are cast in 
the strong conservation form. A fully implicit, finite- 
difference method is used to advance the solution of 
the governing equations in time. A Newton- Rap hson 
subiteration scheme is used to reduce the lineariza- 
tion and factorization errors at each time step. The 


convective terms are evaluated using a third-order- 
accurate upwind-biased Roe scheme. The viscous 
terms are evaluated using second-order accurate cen- 
tral differences and the scheme is second-order accu- 
rate in time. Details of the solution procedure and 
boundary conditions are discussed in Ref. [20] 

TURBULENCE MODELS 

Two models were used to simulate the effects of 
turbulence. The first model is a two-layer algebraic 
turbulence model based on the work of Baldwin and 
Lomax (BL) [21]. Several modifications were made to 
the original BL model based on previous experiences 
with compressor and turbine geometries: 

• The switchover location between the inner and 
outer models cannot move more than a specified 
number of grid points between adjacent stream- 
wise locations. This eliminates non-physical gra- 
dients in the turbulent viscosity near separation 
points. 

• A second derivative smoothing function is used 
on the turbulent viscosity field in separated flow 
regions. This also helps remove non-physical gra- 
dients in the turbulent viscosity in separation re- 
gions. 

• A cutoff value is imposed on the turbulent viscos- 
ity (nominally 1200 times the free stream lami- 
nar viscosity). 

The second model is a two-equation k — e turbu- 
lence model based upon the work of Towne et al. [22]. 
In the current implementation, the equations for the 
turbulence kinetic energy and dissipation rate are 
decoupled from the flow equations and solved using 
an alternating-direction implicit integration scheme. 
Since Newton- Raphson subiterations are used at each 
global time step of the flow solver, decoupling the 
k — e equations from the flow solver should not affect 
the time accuracy of the analysis. The convective 
fluxes in the turbulence equations were discretized us- 
ing first-order accurate upwind differences, while the 
dissipation terms were discretized using second-order 
accurate central differences. The k — e subroutines 
were constructed in a modular manner to allow the 
use of different low Reynolds number approximations; 
the Chien low Reynolds number model has been used 
in the current investigation [23]. 

Note, the k — e equations are solved on both the 
0- and H-grids, while the B-L model is solved only 
on the O-grids. 
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TRANSITION MODEL 


The low Reynolds number environment in low- 
pressure turbines suggests that the flow may be tran- 
sitional. In the current investigation, natural transi- 
tion is modeled using the Abu-Ghannam and Shaw 
(ABS) model [24]. The ABS model determines the 
start of transition based on the momentum thickness 
Reynolds number: 


Re e 


163 4- exp 


F(X e ) 


w; 

6.91 


(1) 


where 


F( A,) = 
F(\,) = 

\ff = 


6.9l4-12.75A e +63.64(A 0 ) 2 A* < 0(2) 
6.91 + 2.48A0 — 12.27 (A^) 2 A* > 0 (3) 


9 2 dUoo 
v dS 


(4) 


This model is considered valid up to turbulence lev- 
els of r — 10%. The end of the transition region is 
calculated as 


points. Again, the average value of y + was 0.20. The 
boundary layers (on average) were discretized with 
25 grid points (see Fig. 2). The dimensions of the O- 
grid was arrived at based on the value of and the 
number of points within the boundary layer, while 
the dimensions of the H-grid was determined by per- 
forming wake convection simulations in the absence 
of airfoils. The finer grid yielded better wake resolu- 
tion; therefore, unless otherwise specified the results 
in the following sections will be based on the finer 
grid. 

The turbine stage was constructed using the same 
airfoil section as was used in the cascade simulations. 
A l-vane/2-blade airfoil count ratio was selected for 
the simulations to reduce the propagation of pres- 
sure waves. The computational grid topology for each 
blade passage was similar to that for the cascade and 
contained 75,393 grid points. The average value of 
y+ was approximately 0.50 for the vane and blade 
airfoils. 

RESULTS 


Re L = 31.8(J2e«) 1,6 ( 5 ) 

In the region between the start and end of transition 
the intermittency function, a, is determined using the 
model developed by Dhawan and Narasimha [25], 

a — 1 — exp (— 4.64£ 2 ) (6) 

Note, the intermittency function, which has a value of 
a = 0 for laminar flow and cr — 1 for turbulent flow, 
is multiplied by value of the viscosity calculated in 
the turbulence model. For cases involving separation 
bubbles the model developed by Roberts [26], and 
modified by Davis et al. [27], is used: 

Re e = 250001o glo (1/ tanh(0.173205r)) (7) 

QEQMETRY AND GRID 

Two geometries have been studied in this investi- 
gation; one is a low-pressure turbine cascade, while 
the other is a low-pressure turbine stage. 

The turbine cascade airfoil design is typical of the 
rotors found in modern aircraft engines (see Fig. 1). 
Two different grid densities were used to discretize 
the turbine cascade. In the coarser discretization the 
O-grid contained 251 x 51 (streamwisex tangential) 
grid points and the H-grid contained 180 x 45 grid 
points, for a total of 20,901 points. The average value 
of y + , the non-dimensional distance of the first grid 
point above the surface, was 0.2. In the finer dis- 
cretization the O-grid contained 281 x 51 and the H- 
grid contained 240 x 45 grid points, for a total 25,131 


Turbine Cascade 

The Mach number at the cascade inlet is M\ — 
0.0897, the inlet flow angle is 35° (measured from 
the axial direction) and the pressure ratio across the 
cascade is P 2 / Pti — 0.9844. The pitch-to-chord ratio 
of the cascade is 0.8856. The Reynolds number, based 
on the inlet velocity and the axial chord of the airfoil, 
was varied from Re — 40, 000 to Re — 120, 000. 

Figures 3 to 8 illustrate the instantaneous Mach 
number and entropy contours predicted in the lam- 
inar, turbulent and transitional simulations, respec- 
tively, at Re = 40, 000 and using the Baldwin- Lomax 
turbulence model. The laminar solution did not 
reach a steady or time-periodic state. The boundary 
layer on the suction side of the airfoil was observed 
to intermittently shed large vortical structures (see 
Figs. 3 and 6). The turbulent solution reached a time- 
periodic state, but still displayed a significant amount 
of unsteadiness (see Figs. 4 and 7). The laminar and 
turbulent entropy contours (see Figs. 6 and 7) also 
show small variations in the entropy upstream of the 
blade passage. The transitional solution reached a 
time-periodic state with much less unsteadiness (see 
Figs. 5 and 8). Two interesting phenomena were ob- 
served in the transitional simulations: 

• Flow transition occurred within the suction sur- 
face separation bubble (near the trailing edge) 
at Reynolds numbers up to Re « 100,000. 

• The location of transition varied in a limit cycle 
(with time) between the beginning of the sepa- 
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ration bubble and aft portion of the separation 
bubble. 

Figures 9 and 10 display comparisons of the isen- 
tropic Mach number distributions predicted in the 
turbulent and transitional simulations, respectively, 
with the design intent for the blade. The results from 
both the coarse and fine grid topologies are included. 
Both the turbulent and transitional solutions exhibit 
good agreement with the design data, although the 
turbulent solution shows more variation between the 
coarse and fine grids. Table 1 contains time-averaged 
flow quantities at the inlet and exit of the cascade. 
The losses, AP*/P<i, are lowest in the laminar flow 
simulation and greatest in the turbulent simulation. 
The average inlet/exit flow angles and Mach numbers 
exhibit good agreement with the design intent. 

Simulations were also performed using the k — e 
turbulence model. In the simulations the free stream 
turbulence intensity wa is varied between r — 3% and 
r — 6%. In lieu of a dissipation length scale, the cur- 
rent implementation of the k — e turbulence model 
requires a value for the free stream turbulent viscos- 
ity (from which the length scale can be obtained). 
Unless specified, the value of free stream turbulent 
viscosity was set to 10 times the free stream laminar 
viscosity. Figure 1 1 illustrates instantaneous entropy 
contours at Re = 80, 000. The flow field achieves a 
nearly steady-state condition at this Reynolds num- 
ber. A comparison of the predicted and design isen- 
tropic Mach number distributions is shown in Fig. 12. 
There is good agreement between the two distribu- 
tions. Figures 13 shows the variation of the skin fric- 
tion as the free stream turbulence intensity level is 
varied. The solutions are similar, except near the 
leading edge on the suction surface where the skin 
friction decreases with increasing turbulence inten- 
sity. Figure 14 shows the variation of the skin fric- 
tion as the free stream turbulent viscosity (dissipa- 
tion length scale) is varied. It is observed that the 
solution is very dependent on the length scale. This 
underscores the need for experimental guidance in de- 
termining the dissipation length scale. 

Several wake-passing simulations were performed 
using a hyperbolic tangent distribution to simulate 
the wake profile. The wake was moved across the 
inlet of the cascade using the boundary conditions 
outlined by Dorney and Verdon [20]. There was one 
wake specified for two blade passages and the velocity 
deficit was set at 5% of the average inlet velocity. The 
Reynolds number was specified to be Re — 120,000 
and the flow was assumed to be transitional. Fig- 
ures 15-22 show instantaneous entropy contours and 
momentum thickness distributions at 0%, 25%, 50% 
and 75% of a wake-passing cycle, respectively. The 


momentum thickness distributions are for the lower 
blade in the cascade. At the beginning of the cycle 
(see Figs. 15 and 16) the wake has just passed over 
the leading edge and has started moving along the 
pressure surface of the blade. The corresponding mo- 
mentum thickness distribution displays elevated val- 
ues from approximately 20% to 60% of the axial chord 
on the pressure surface. The relatively large values of 
the momentum thickness on the pressure surface are 
caused by the existence of a cove separation bubble. 
As the wake is convected along the pressure surface 
of the blade (see Figs. 17 and 18), the values of the 
momentum thickness are decreased in the mid-chord 
region. At 50% of the wake-passing cycle (see Figs. 19 
and 20), the wake is being slowly convected along the 
pressure side of the blade passage. The wake is be- 
ing stretched from the pressure surface of one blade 
to the suction surface of the adjacent blade because 
of differences in the convection speed. During this 
time the peak value of the momentum thickness con- 
tinues to decrease, and elevated values are observed 
over only 30% of the pressure surface. Thus, the 
size of the cove separation bubble is decreased by 
the passing of the wakes. At 75% of the cycle the 
next wake has just impinged on the leading edge of 
the blade (see Figs. 21 and 22). The corresponding 
momentum thickness distribution again displays ele- 
vated values on the pressure surface near mid-chord. 
Throughout the wake-passing cycle the suction sur- 
face values of the momentum thickness varied only a 
small amount. The small variations were caused by 
unsteady changes in the location of (natural) transi- 
tion. 

Turbine Stage 

Numerical simulations were performed for the stage 
configuration (see Fig. 23) at Reynolds numbers from 
Re = 27, 000 to Re = 120, 000 (based on upstream 
flow quantities). The inlet flow conditions for the 
stage were similar to those for the cascade. The ro- 
tational speed of the rotor was chosen such that the 
time- averaged relative flow angle entering the rotor 
passage was approximately equal to that used in the 
isolated blade row simulations (35°). Both turbulent 
(BL model) and transitional flow simulations were 
conducted on the stage geometry. Two techniques 
were used to perform the transitional simulations: 

• The location of transition was allowed to “float” 
as the solution advanced in time (i.e., the tran- 
sition model is used to calculate the location of 
transition at each time step). 

• The location of transition was “fixed” at the 
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time-averaged position determined from the 
floating transition simulation. 

Figure 24 illustrates the variation in the location of 
transition on the suction surface of the vane and ro- 
tor as function of the Reynolds number. In the stage 
simulations the transition from laminar to turbulent 
flow always occurred as bubble transition, and usu- 
ally occurred near the trailing edge. 

Figures 25 to 28 illustrate instantaneous Mach 
number and entropy contours at Reynolds numbers 
of Re = 40, 000 and Re = 120, 000, respectively, for 
turbulent flow conditions. The Mach and entropy 
contours at Re = 40,000 indicate thicker boundary 
layers and more unsteadiness than at Re = 120, 000. 
The entropy contours at Re = 40, 000 also indicate 
the presence of waves upstream of the vane passage. 
Figures 29 and 30 illustrates the unsteady pressure 
envelopes at the two Reynolds numbers. Under tur- 
bulent flow assumptions the two sets of pressure en- 
velopes are similar with the exception that there are 
somewhat larger excursions about the time-averaged 
values at Re = 40, 000. 

Figures 31 to 36 display skin friction envelopes un- 
der the turbulent, floating transition and fixed transi- 
tion flow assumptions, respectively, for Re = 40, 000 
and Re = 120,000. Assuming turbulent conditions 
(see Figs. 31 and 32) the flow remains attached on 
the vane at both Reynolds numbers. There is a cove 
separation bubble on the pressure surface of the rotor 
at both Reynolds numbers, but the extent of the sep- 
aration is much greater at Re = 40,000. Under the 
turbulent flow conditions, the efficiency of the stage 
at Re = 120, 000 is approximately 1 point higher that 
at Re — 40,000 (see Tables 2 and 3). Note that due 
to the low Mach number of the flow small changes in 
the total temperature and/or total pressure can result 
in a large change in the efficiency. In the case of the 
floating transition (see Figs. 33 and 34), the solution 
at Re = 40,000 shows significant unsteadiness and 
flow separation. The separated flow regions include: 
1) the aft 25% of the vane suction surface, 2) the aft 
25% of the rotor suction surface and 3) the first 40% 
of the rotor pressure surface (cove separation). These 
large separated flow regions result in a very low stage 
efficiency (see Table 2). At Re = 120,000 the float- 
ing transition simulation displays much less unsteadi- 
ness. The flow separates on the suction surface of 
the vane, but quickly reattaches. The time-averaged 
solution on the rotor suction surface contains only 
a small separated flow region, although the size of 
the separation region alternately grows much larger 
and disappears completely during a blade-passing cy- 
cle. Using the floating transition assumptions the ef- 
ficiency at Re = 120, 000 is 10 points higher than 


at Re = 40, 000 due to much smaller separated flow 
regions (see Tables 2 and 3). In the case of fixed tran- 
sition, the extent of the flow separation on the suc- 
tion surface of the vane and blade are substantially 
reduced at both Reynolds numbers. The cove separa- 
tion on the pressure surface of the rotor, however, is 
slightly larger at both Reynolds numbers. In general, 
there is much less unsteadiness when the transition 
locations are fixed at their time-averaged positions. 
The stage efficiency at Re = 120, 000 is still 10 points 
higher than at Re = 40,000, and, surprisingly, at 
both Reynolds numbers the efficiency is lower when 
the location of transition is fixed than when it is al- 
lowed to float (see Tables 2 and 3). 

Figures 37 to 42 display the momentum thickness 
envelopes under the turbulent, floating transition and 
fixed transition flow assumptions, respectively, for 
Re ~ 40,000 and Re = 120,000. The momentum 
thickness distributions, in general, follow the same 
trends as the skin friction. The interesting feature 
of the Figs. 37 to 42 is the size of momentum thick- 
ness in the cove region on the pressure surface of the 
airfoils. The momentum layer is thicker and exhibits 
more unsteadiness in this region than over the ma- 
jority of the suction surface. 

CONCLUSIONS 

A series of cascade and stage simulations have been 
performed for a low-pressure turbine geometry. The 
numerical models used to simulate the physical mech- 
anisms, as well as the physics associated with low 
Reynolds number flows, have been investigated. Both 
the cascade and stage simulations revealed that the 
predicted losses and efficiency are highly dependent 
on the assumed state of the boundary layers ( i.e 
laminar, transitional or turbulent). In the transi- 
tional simulations the location of transition was in- 
fluenced by both the periodic (wake or blade passing) 
and random (boundary layer shedding) unsteadiness 
in the flow. Work is in progress to incorporate more 
sophisticated transition models into the flow analysis. 
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Turb. 

Float Trans. 

Fix Trans. 

M 1 

0.0897 

0.0848 

0.0844 

m 2 

0.1409 

0.1357 

0.1338 

m 3 

0.0914 

0.0858 

0.0840 

m 2 

0.0909 

0.0859 

0.0839 

m 3 

0.1422 

0.1366 

0.1344 

<*1 

—35.00° 

-35.00° 

-35.00° 

«2 

57.71° 

58.17° 

58.27° 

a 3 

-34.67° 

-33.95° 

-32.84° 

ft 

33.29° 

32.33° 

32.20° 

ft 

-58.15° 

-58.74° 

-58.43° 

Vtt 

0.925 

0.847 

0.826 


Table 2: Time-averaged inlet and exit flow quantities 
- stage - Re=40,000. 



Turb. 

Float Trans. 

Fix Trans. 

Mi 

0.0903 

0.0902 

0.0890 

m 2 

0.1419 

0.1431 

0.1410 

M 3 

0.0919 

0.0923 

0.0907 

m 2 

0.0916 

0.0923 

0.0903 

Ms 

0.1429 

0.1437 

0.1419 


-35.00° 

-35.00° 

-35.00° 

<*2 

57.82° 

58.24° 

58.25° 

a 3 

-34.88° 

-35.53° 


ft 

33.88° 

34.90° 

34.40° 

ft 

-58.17° 

-58.49° 

-58.51° 

Vtt 

0.934 

0.947 

0.929 



Design 

Lam. 

Turb. 

Trans. 

Mi 

0.0897 

0.0907 

0.0927 

0.0934 

m 2 

0.1500 

0.1474 

0.1440 

0.1468 

<*i 

35.00° 

35.00° 

35.00° 

35.00° 

<*2 

-60.00° 

-58.65° 

-58.65° 

-58.27° 

AP,/P t i 

— 

0.00122 

0.00063 

0.00082 


Table 1: Time-averaged inlet and exit flow quantities Table 3: Time-averaged inlet and exit flow quantities 
- cascade - Re=40,000. - sta S e " Re=120,000. 
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Figure 1: Low-pressure turbine cascade. 


Figure 3: Instantaneous Mach contours - Re=40,000 
- laminar. 



Figure 2: Velocity profiles on blade suction surface - j?ig Ure 4 ; Instantaneous Mach contours - Re=40,000 
Re 40,000 - laminar. _ turbulent - BL. 








Figure 5: Instantaneous Mach contours - Re=40,000 Figure 7: Instantaneous entropy contours 

- transitional - BL. Re=40,000 - turbulent - BL. 



Figure 6: Instantaneous entropy contours - Figure 8: Instantaneous entropy contours - 

Re=40,000 - laminar. Re=40,000 - transitional - BL. 
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Figure 9: Isentropic Mach number distribution - Figure 11: Instantaneous entropy contours - 

Re=40,000 - turbulent - BL. Re=80,000 - k — e - Tu — 3% - a*too = 10.0. 
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Figure 10: Isentropic Mach number distribution - Figure 12: Isentropic Mach number distribution 
Re=40,000 - transitional - BL. Re=80,000 - k — e- Tu — 3% - utoo — 10.0. 
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Figure 17: Instantaneous entropy contours at 25% Figure 19: Instantaneous entropy contours at 50% 
wake-passing cycle - Re=120,000. wake-passing cycle - Re=120,000. 



Figure 18: Momentum thickness distribution at 25% Figure 20: Momentum thickness distribution at 50% 
wake-passing cycle - Re=120,000. wake-passing cycle - Re=120,000. 






Figure 21: Instantaneous entropy contours at 75% 
wake-passing cycle - Re=120,000. 



Figure 22: Momentum thickness distribution at 75% 
wake-passing cycle - Re= 120,000. 


Figure 23: Low-pressure turbine stage. 



Figure 24: Time-averaged transition locations. 





Figure 25: Instantaneous Mach contours - Re=40, 000 Figure 27: Instantaneous entropy contours - 

- turbulent - BL. Re=40,000 - turbulent - BL. 



(S/Soo) roin =0.999 
(5/Sco) max = 1.009 



Figure 26: Instantaneous 

Re=120,000 - turbulent - BL. 


Mach contours - Figure 28: 
Re=120,000 


Instantaneous entropy contours 
turbulent - BL. 
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Figure 29: Unsteady pressure envelope - Re=40,000 Figure 31: Skin friction envelope - Re— 40,000 - tur- 
- turbulent. bulent. 



Figure 30: Unsteady pressure envelope - Re=120,000 
- turbulent. 


Figure 32: Skin friction envelope - Re=120,000 - tur- 
bulent. 
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Figure 33: Skin friction envelope - Re=40,000 - float- 
ing transition. 



Figure 35: Skin friction envelope - Re=40,000 - fixed 
transition. 



Figure 34: Skin friction envelope - Re=120,000 - 
floating transition. 



Figure 36: Skin friction envelope - Re=120,000 - fixed 
transition. 
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Figure 37: Momentum thickness envelope 

Re=40,000 - turbulent 


Figure 39: Momentum thickness envelope - 

Re=40,000 - floating transition 




Figure 38: Momentum thickness envelope 

Re=120,000 - turbulent 


Figure 40: Momentum thickness envelope - 

Re= 120,000 - floating transition 
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